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We study the phase structures of charge neutral quark matter under the /3-equilibrium for a wide 
range of the quark-quark coupling strength within a four-Fermion model. A comprehensive and 
unified picture for the phase transitions from weak to strong coupling is presented. We first develop 
a technique to deal with the gap equation and neutrality constraints without recourse to numerical 
derivatives, and show that the off-diagonal color densities automatically vanish with the standard 
assumption for the diquark condensates. The following are shown by the numerical analyses: (i) 
The thermally-robustest pairing phase is the two-flavor pairing (2SC) in any coupling case, while 
the second one for relatively low density is the up-quark pairing (uSC) phase or the color-flavor 
locked (CFL) phase depending on the coupling strength and the value of strange quark mass, (ii) If 
the diquark coupling strength is large enough, the phase diagram is much simplified and is free from 
the instability problems associated with imaginary Meissner masses in the gapless phases, (iii) The 
interplay between the chiral and diquark dynamics may bring a non-trivial first order transition even 
in the pairing phases at high density. We confirm (i) also by using the Ginzburg-Landau analysis 
expanding the pair susceptibilities up to quartic order in the strange quark mass. We obtain the 
analytic expression for the doubly critical density where the two lines for the second order phase 
transitions merge, and below which the down-quark pairing (dSC) phase is taken over by the uSC 
phase. Also we study how the phase transitions from fully gapped states to partially ungapped 
states are smeared at finite temperature by introducing the order parameters for these transitions. 

PACS numbers: 12.38.-t, 25.75.Nq 



I. INTRODUCTION 

Surprisingly rich phase structure of quark matter is 
being revealed by extensive theoretical studies [1-3]. On 
the basis of the asymptotic-free nature of QCD and the 
attraction between quarks due to a gluon exchange, it 
is now believed that the ground state of the sufficiently 
cold and extremely dense matter is in the color-flavor 
locked (CFL) phase where all the quark species equally 
participate in pairing [4]. In contrast, it remains still 
controversial what phase sets in next to the CFL as the 
density is decreased and/or the temperature is raised. 

At zero temperature (T = 0), the following two in- 
gredients play crucial roles for determining the second 
densest phase in QCD; (i) the strange quark mass [5- 
7], and (ii) the charge neutrality constraints as well as 
the /3-cquilibrium condition [8-10]. The former tends to 
bring a Fermi-momentum mismatch between the light 
and strange quarks, while the latter in the paired phase 
tends to match the Fermi-momenta by tuning the charge 
chemical potentials. Some theoretical studies suggest 
that, when the baryon density is decreased, the effect 
(ii) in the CFL phase cancels the effect (i) down to some 
critical density, at which the CFL phase turns into the 
gapless CFL (gCFL) phase [11]. In the realistic quark 
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matter in compact stars, the effect of (ii) gives so strong 
constraint because of a long range nature of gauge inter- 
actions that such an exotic phase may exist stably. The 
gCFL phase has been extensively studied and claimed to 
be the most promising candidate of next phase down in 
density at vanishing temperature [11] and is also known 
to lead an interesting astrophysical consequence on the 
cooling history of an aged compact star [12]. It should 
be, however, noted here that the gapless phases are usu- 
ally accompanied by some transverse gluonic modes with 
imaginary Meissner mass at low temperature [13, 14], in- 
dicating the existence of more stable exotic states [15-21]. 
Although the resolution of the instability or finding the 
possible new stable state is now one of the central issues 
in this field, we won't deal with this problem; neverthe- 
less we will give some suggestions to possible resolution 
on the basis of the results obtained in this work. 

The neutrality constraints are also known to bring 
about an interesting phase even at finite temperature 
(T ^ 0). In fact, a Ginzburg-Landau analysis shows that 
the d quark pairing phase with the u-d and d-s pairings, 
denoted by "dSC" , appears as the second phase when the 
2SC (CFL) phase is cooled (heated) [22]. The existence 
of such a kind of intermediate phase has also been con- 
firmed using the NJL model [23] . Also the gapless version 
of the dSC with an additional gapless mode is also exam- 
ined for finite temperature [22] and for zero temperature 
[24]; the gapless dSC is shown to be always a metastable 
at zero temperature, while it may exist stably for nonzero 
temperature. It should be noted, however, that the gap- 
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less phase for nonzero temperature cannot be thcrmo- 
dynamically distinguished from the fully gapped one; in 
fact, it has no distinct phase boundary with the gapped 
one, and thus is of less physical interest than that at zero 
temperature. 

In the most of the literature concerning the pairing 
phases of the quark matter [11, 23, 25], the strange 
quark mass is treated as a parameter just like an ex- 
ternal magnetic field applied to a metallic superconduc- 
tor [26]. In QCD, however, the strong attraction exists 
in the scaler quark-antiquark channel which leads to a 
non-pcrturbativc phenomenon called the dynamical chi- 
ral condensation in the low density regime. It is inter- 
esting to investigate how the incorporation of the dy- 
namical formation of the chiral condensates affects color 
superconducting phases [27, 28]. The competition be- 
tween the chiral and diquark condensations under the 
charge neutrality constraints was first investigated in a 
four-Fermion model [24]; it was shown that the second 
densest phase next to the CFL strongly depends on the 
quark-quark (diquark) coupling strength, and in partic- 
ular, the gapless CFL phase may be washed out from the 
phase diagram except for an extremely weak coupling 
case. Their investigation is extended to finite tempera- 
tures [29, 30] and recently also applied to the system with 
finite neutrino density [31]. Some of the results obtained 
in [29, 30], however, seemingly contradict with the pre- 
vious claim by the authors of [22, 23] on the appearance 
of the dSC phase. This is one of puzzles indicating that 
a further investigation is needed for a systematic under- 
standing of the phases of quark matter and transitions 
among them. 

A main aim of this paper is to explore the phase di- 
agram for a wide range of diquark coupling systemat- 
ically and thereby give a unified picture for the phase 
transitions from weak to strong coupling. We clarify the 
detailed mechanism and features of thermal phase tran- 
sitions, and make clear the relations among the previous 
studies [22, 23, 29, 30]. In particular, we shall study 
the the following, (i) How the (/x, T)-domain accommo- 
dating the pairing grows and that for the gapless phases 
vanishes as the coupling in the scaler diquark channel is 
increased, (ii) How we can draw a unified picture for the 
thermal unpairing phase transitions from weak to strong 
coupling; for this purpose, we extend the previous study 
[22] by making re-analysis of Ginzburg-Landau expan- 
sion focusing on the effects of the terms quartic in the 
strange quark mass, (iii) How and why non-trivial first 
order phase transitions could be caused by the compe- 
tition between the pairing and chiral dynamics in the 
strange quark sector, (iv) How the transition from the 
fully gapped phase to the partially ungapped (gapless) 
one is smeared at finite temperature by introducing a 
definite order parameter. 

The paper is organized as follows. In Sec. II, we in- 
troduce the model and derive the gap equation under 
the neutrality constraints. Some technical developments 
are presented, by which the numerical derivatives can be 



replaced by simple algebraic equations. In Sec. Ill, we 
discuss the points listed above based on the numerical 
results. Also the Ginzburg-Landau analysis is performed 
and some interesting aspects of the thermal phase transi- 
tions are given at the end of this section. The summary 
and outlook are provided in Sec. IV. In Appendix A, 
we derive the Ginzburg-Landau potential expanded up 
to quartic order in the strange quark mass. In Appendix 
B, we show that the off-diagonal color densities automat- 
ically vanish under the standard ansatz for the diquark 
condensates. 



II. FORMULATION 

In this section, we introduce our model and formu- 
late the mean field approximation under the kinetic con- 
straints. In Sec. II B, we present some useful algebraic 
techniques to solve the gap equation and constraints 
without recourse to numerical derivatives. 



A. Model 

We start with the following Lagrangian as in [24] 

3 

r?=l 

+ ^-[(q\ F q) 2 + &J 5 \ F q) 2 }. (1) 

Here, = {-\/2/31, \f} are the unit matrix and the 
Gell-Mann matrices in the flavor space. P v is defined as 
[11] 

(P v )1j = ij 5 Ce vab e vl j no sum over index r\ (2) 

and P v = 7o(Pt,) 7o- a , b, ■ ■ ■ and ■ ■ ■ represent the 
color and flavor indices, respectively. The second term 
in Eq. (1) simulates the attractive interaction in the 
color anti-triplet, flavor anti-triplet and J p = + channel 
in QCD. rag = diag.{7« u , m<2, m s } is the current-quark 
mass matrix. In order to impose the color and electric 
neutralities, we have introduced in Eq. (1) the chemical 
potential matrix fi in the color-flavor space as 

(t$ = fi - fieQa + ^n ib + vsTt 

+(XiT? b + inTf + ^Tt (3) 

Q = diag.(2/3, — 1/3, — 1/3) counts electric charge of 
each quark species. T 3 = diag.(l/2, — 1/2, 0) and T 8 = 
diag.(l/3, 1/3, —2/3) are the diagonal charges for quarks 
in the fundamental representation of the color SU(3). We 
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have included the second and third lines in Eq. (3) which 
stand for the chemical potentials for off-diagonal color 
charges: It has been recently claimed that these chemical 
potentials should be included for the complete color neu- 
tralities [32]. We can prove, however, that if the diquark 
condensates have the color-flavor structure as given by 
Eq. (2), then the off-diagonal color densities must auto- 
matically vanish for Mi = M2 = M4 = M5 = M6 = M7 = 0; 
see Appendix B for the detail of the proof. Thus, we 
can safely adopt the usual diagonal ansatz for the chem- 
ical potential matrix. The explicit forms of the diagonal 
elements (Mai = (if?) are as follows. 
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We treat the diquark coupling constant Gd as a sim- 
ple parameter, although the perturbative one-gluon ex- 
change vertex £ int = -(g 2 /2)qj ll (\ a /2)qq^(\ a /2)q, 
which is valid at extremely high density, tells us that 
Gd/G s = 1/2 with N c = 3 [33-35]. As a measure of cou- 
pling constant Gd, we shall mainly use the gap energy 
(Ao) in the pure CFL phase at (1 = 500 MeV and T = 
in the chiral SU(3) limit, as in [11, 23, 24]. 

We evaluate the thermodynamic potential in the mean- 
field approximation; 
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are the gap parameter and constituent quark mass ma- 
trices. 5 denotes the 72 x 72 Nambu-Gor'kov propagator 



defined by 
S' _1 (iuj n ,p) 



(8) 



with j> = iuj n "fo ~ P ' 1- Finally, fi e is the contribution 
from massless electrons 
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The optimal values of the variational parameters A n , M 
and M s must satisfy the stationary condition (the gap 
equations); 
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Our task is to search the minimum of the effective po- 
tential by solving these gap equations under the local 
electric and color charge neutrality conditions; 

<9fi dfi , 9fi 

^- = 0, — = and — = 0. (11) 

For a later convenience, we define here the electron den- 
sity by 
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The formulation made above is a straightforward exten- 
sion of our previous work [24] to the (T ^ 0) case [29, 30]. 



B. Representation of gap equation and kinetic 
constraints in terms of quasi-quark wave functions 

In this section, we present some analytical way to deal 
with the gap equation and neutrality constraints. We 
shall show that the gap equation and neutrality con- 
straints derived above can be further simplified with the 
aid of the quasi-quark wave functions (spinors). By doing 
this, not only the physical meanings of these equations 
become transparent, but also the numerical derivatives 
can be circumvented. In particular, the latter has a prac- 
tical advantage because the computations of matrix ele- 
ments arc more favorable than the numerical derivatives. 

First, we introduce the Nambu-Gor'kov mean field 
Hamiltonian density TC(p; (1, A, M) following [11, 23], 



T S (iu) n ,p)T ■■ 
Here we have defined 



iu n l 72 -H(p;(i,A,M). (13) 
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i 36 
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I72 denotes the 72-dimensional unit matrix and TL is 
Hamiltonian density which is also 72 x 72 matrix in 
the color, flavor and spinor space. In the same man- 
ner as shown in [29], we can lift the spin degeneracy 
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away from Tt as H = Ji^Q (p) (g) P + © Ti^Q (— p) <g> P_ give the explicit form of these matrices below, 
with P± = 1± P' rT being the helicity projectors. H-^ip) 



has a block-diagnalized form like H.3e(p) = W12 
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Gor'kov degeneracy is manifest in the latter three blocks; 
H ( s a ' 0) = U { t P) ®[-n { t P) ] ■ Thus we need to know only 
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Ti.i2 S f° r the evaluation of the effective potential. We 



Uf and H[ u 2 r ' d9 > sb) 



(M a -iJL a p -iA a P \ 

p -M a - p a -iA a(3 
iA af} Mp + pp -p 

V -p -Mp + pp J 

(15) 

where A a/3 = Ai, A 2 and A3 for (a, (3) = (bd, sg), (sr, ub) 
and (ug,dr), respectively, and 
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Because the Nambu-Gor'kov degeneracy is not removed, 
the eigenvalues of this matrix have six sets of the doublet 
(e, — e) corresponding to the energies for a quasi-quark 
and its Nambu-Gor'kov partner (anti-quasi-quark). 

Let us write 18(= 72 -j- 2 -7- 2) independent eigen- 
values of H as {e Q } with the index a redefined to 
the label for the quasi-quarks (a = 1,2, ••• ,18) from 
that for the direct products of color and flavor (a = 
ur, dg, sb, db, sg, sr, ub, ug, dr). Then the thermodynamic 
potential can be simplified to 
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We have subtracted the vacuum contribution in the sys- 
tem with the nine massless quarks. It is difficult to nu- 
merically search the minima with respecting the neutral- 
ity constraints from this effective potential alone. There- 
fore, we search them with the aid of the gap equation 
which is the A-derivative of the effective potential. If 
possible, numerical derivatives should be avoided because 
the numerical errors associated with them are not well 
controllable. To avoid them, we express the gap equa- 
tions in terms of the eigenvectors (cigen-spinors) of Ti.. 
First, we simply differentiate Eq. (5) with respect to A^ 



and equating the result to zero to obtain, 
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Here H36 is the reduced Nambu-Gor'kov Hamiltonian 
density defined by removing the spin degeneracy from H. 
as was introduced above. Also we note that the reduced 
Hamiltonian density takes the form 



n 36 = Ho + M u <f> Mu + M d <f> Md +M s cj) 
+Ai0 x + A 2 4> 2 + A 3 3 
-3fiB - p e Q e - fx 3 T 3 - /z 8 T 8 



(18) 



with Hq = 7i36|^=o,A=o,Af=o being the Hamiltonian 
density for the system with nine free massless quarks. 
4>Mi i 't'Ar, an d B, Q e , T3 : 8 are the 36 x 36 matrices with 
constant matrix elements. These matrices can be ob- 
tained by differentiating the reduced Hamiltonian matrix 
with respect to M u> d, M s , Ai^, ^, Me, 3. 8 as follows. 
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Because 7i 3 6 is an Hermitian matrix for any momen- 
tum, we can define the complete set of spinors {\p, a, ±)} 
by the eigenvalue equation 



H 36 \p,ot,±) = ±e a {p)\p,a,±) 



(20) 
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a = 1, • ■ ■ ,18 and ± denote the Nambu-Gor'kov charges. 
Then we have 

de 

(p,a,±\<t> n \p,a,±) = ±-^-, ( V = 1,2,3). (21) 
Thus we find that Eq. (17) is reduced to 

18 

a—1 K ' 



where the Matsubara summation has been performed. In 
much the same way, we obtain 
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These formulae enable us to evaluate the gap equation 
solely by computing the 18-dimensional eigen-spinors de- 
fined by Eq. (20) and some matrix elements in these bases 
[36] without recourse to numerical derivatives done in 
[11, 23, 29]. 

The charge neutrality constraints can be also expressed 
in terms of the matrix elements in the basis composed 
of these eigen-spinors. A straightforward application of 
the above method to p 3>8 = —dVl q /dpz$ = 0, p\ = 
—dil q /dp e = —n e and pb = — \dVl q / dp leads 
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Here, the Fermi-Dirac distribution function fp (e) = 
l/(e e / T + 1) is introduced. 

At zero temperature, /f(e) = 0(— e) so that the net 
charge will be accumulated in the blocking region where 
e < 0. For the neutrality constraints to be satisfied, there 
must be the opposite charge density coming from a non- 
quark sector or from the quark sector with the finite back- 
ground charge density / £a(P; a > +\Q\p, a , +> ^ 0, 
which are supplied by tuning the charge chemical poten- 
tials p e , pz and pg. 

For a later convenience, we introduce here a charge 
generator 



Q 



-Q-n 
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the operation of which keeps the CFL state invariant 
(neutral) and thus represents an unbroken U(l) symme- 
try in the CFL phase [3]. If we choose the orthogonal 



basis of the broken charges (X, Y) in addition to unbro- 
ken Q as [10], 

X = Q + T 3 -4T 8 , 
Y = Q-T 3 , 
the chemical potentials for these charges become 

= Ts(-MQ +M3 - 4/i 8 ), 

Mr = ±{-pQ-p 3 )- 

We find that the matrix representation of Q operator 
in the 32-dimensional color-flavor mixed Nambu-Gor'kov 
base is given by 

«— ^--«.-^-5 r - < 23 > 

Because this commutes with the Nambu-Gor'kov Hamil- 
tonian density 7^36 as it should be, the quasi-particles 
(eigen-spinors) have definite Q~ c h ar g es which can be 
shown to be of integers (+1,0,-1) [3]. It is also to be 
noted here, that thermodynamic potential in the quark 
sector n q does not depend on pq in the fully gapped 

CFL phase so that it is a Q-insulator [10]. Because of 
this, the value of pq in the CFL phase (T = 0) can- 
not be determined by the Q neutrality condition in the 
quark sector; it should be determined completely by the 
vanishing-point of very gentle slope of potential curvature 
coming from the electron sector (pl/3ir 2 = 0). 



III. NUMERICAL RESULTS AND 
DISCUSSIONS 

In this section, we present our numerical results of the 
solution of the gap equations, and provide the phase di- 
agrams in the (fx, T)-plane for several values of diquark 
coupling G d - 

Before that, we fix our model parameters. We take 
the chiral SU(2) limit for the u, d current quark masses 
(m u = m d = 0) and m s = 80McV [24]. These 
values might slightly underestimate the effect of the 
current masses because m Utd (2GcV) = 3-4 MeV and 
m s (2GeV) = 80-100 MeV according to the full lattice 
QCD simulation [37]. Also we restrict the variational 
space by putting M = M u = Md for simplicity. This 
simplification does not matter in the chiral symmetry re- 
stored phase [29, 30]. 

For comparison with the previous work [29], we write 
down the dimensionlcss parameters adopted in this 
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TABLE I: The nonzero gap parameters, some conditions between gaps and chemical potentials, and the gapless quarks in each 
phase. The figure in the parenthesis in the first column represents the number of gapped quasi-quark mode. We have defined 
the (ai-bj) relative chemical potential by fiaibj = (^ai — /J.bj)/2. u dg-sb (1)" means that one of the linear combinations, the 
dg quark and the sb hole remains gapless. The equation for chemical potentials in a bracket must necessarily hold for some 
symmetry or kinetic reason. 



study; 

m Ujd /A = 0, m s /A = 0.1, 
r a 2 

Gs A 2 = = 2.17. 

SN C 

The value of G s is chosen so that the dynamical quark 
mass at n = is 400 MeV for the cutoff A = 800 MeV just 
for comparison with our previous work [24]. Note, how- 
ever, that this coupling is a little larger than the value 
extracted in the NJL model analysis of the meson spec- 
troscopy with instanton induced six-quark coupling [38] , 
i.e., GsA 2 = 1.835 which is adopted in [29]. We shall per- 
form the calculation with following five different values 
of the diquark coupling (Gd = Gd/16): 

1. extremely weak coupling: 

A = 25 MeV <-► G^A 2 = 0.91 (G D /G S = 0.42) 

2. weak coupling: 

A = 80 MeV ^ G D A 2 = 1.37 (G D /G S = 0.63) 

3. intermediate coupling: 

A = 125 MeV ^ G D A 2 = 1.69 (G D /G S = 0.78) 

4. strong coupling: 

A = 160 MeV <-> G D A 2 = 1.92 (G D /G S = 0.88) 

5. extremely strong coupling: 

A = 200 MeV <-> G D A 2 = 2.17 (G D /G S = 1.00) 

The values of Gd/Gs for the intermediate and extremely 
strong coupling cases are similar to those employed in 
[29], i.e., G D /Gs = 3/4 and G D /G S = 1. However, 



the following notice is in order here, (i) Our strange 
quark mass m s /A = 0.1 is about one half of the value 
m s /A = 0.23 adopted in [29]. (ii) We did not included 
the six-quark interaction which effectively increases the 
scalar coupling G5. These two differences make our case 
favorable to the pairing phases rather than the unpaired 
quark matter (UQM) phase or the chiral- symmetry bro- 
ken (xSB) phase. In fact, we will see that the phase 
diagrams for our weak coupling and intermediate cou- 
pling cases seem more or less to correspond to those for 
G D /Gs = 3/4 and G D /G S = 1 in [29], respectively. 

We consider the candidates of phase listed in TABLE 1 
as in [24] in the numerical analyses below. 



A. Phases for extremely weak coupling 

Let us first discuss the extremely weak coupling case 
(Ao = 25 MeV). The phase structure in this case is dis- 
played in Fig. 1. At a first glance, we notice that the 
UQM phase without any symmetry breaking dominates 
the phase diagram pushing the superconducting phases 
to the high density regime. This is because the energy 
gain due to the (ss) condensation in the UQM phase is 
larger than the paring energy under the stress as is clar- 
ified in [24] . There are also several thermal phase transi- 
tions. The thermally-robustest pairing phase is the 2SC 
and the second phase in this case is the uSC as is found 
in [29, 30]. In the following, we shall discuss the fea- 
tures of these phase transitions in detail. We first make 
a close examination on the zero temperature case, and 
then investigate the finite temperature case. 
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FIG. 1: The phase diagram for the extremely weak coupling 
(Ao = 25 MeV). Sold (dashed) lines mean first (second) order 
transitions. Large dot placed on the /^-axis (n = 581.1 MeV) 
represents the transition point at which the CFL phase on the 
high density side continuously turns into the gCFL phase on 
the low density side. 



a. CFL/gCFL and gCFL/UQM transitions atT = 0: 
The phases realized at T = are the CFL, gCFL and 
UQM states. In Fig. 2(a), we show the gaps, masses, 
and chemical potentials along the (g)CFL solution of the 
gap equation. For [i > 558.4 MeV we confirm that the 
condensation energy in the (g)CFL solution is largest 
among those for all the candidates. The excitation gaps 
of the nine quasi-quarks in the CFL phase are shown in 
FIG. 2(b). We can see from the figure that one quasi- 
quark has the largest excitation gap, and other quarks 
have relatively small gaps. The latter eight modes are the 
remnants of the color-flavor octet modes in the pure CFL 
phase at M s = 0. In the CFL phase, the non-vanishing 
/Us and M s are dynamically realized so that the original 
symmetry of color-flavor diagonal SU(3)c+y is explicitly 
broken down to color S\J(2)c+v (color-flavor isospin) in 
the Lagrangian level; because of this, the octet modes 
split into the isospin-singlet mode (like eta) and two set 
of doublet (kaonic) modes and the triplet (pionic) modes 
as 8 -► 3 (+1, 0, -1) + 2 (+1, 0) + 2 (0, -1) + 1 (0) where 
the associated Q-chargc is indicated in the parenthesis. 
One of the most striking features of the CFL phase is 
the absence of electrons; the electric neutrality is real- 
ized solely by the quark sector n u = = n s as can be 
seen in Fig. 2(c). This CFL phase behaves as Q-insulator 
because of the absence of gapless Q-carriers [11]. 

As the density is decreased, the stress energy 
—Hs(iJ,)/2 + M^{ji)/Afi [11] becomes large in the CFL 
state. When it reaches Ai = A2, the first qualitative 
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FIG. 2: (a) The (g)CFL solution of the gap equation under 
the neutrality constraints. The (g)CFL is energetically taken 
over by the UQM phase on the left side of \i = 558.4 MeV 
denoted by the vertical dashed line, (b) The excitation gaps 
in the nine quasi-quark spectra as a function of /i. The Q- 
charge associated with each mode is also indicated, (c) The 
quark and electron densities are displayed. The lines on the 
left hand side are densities in the UQM phase, while those on 
the right hand side are densities along the (g)CFL solution. 
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change takes place; this happens when the chemical po- 
tential is decreased down to /i = 581.1 MeV. At this 
point, the CFL state continuously turns into the gCFL 
phase. Just at the gCFL onset, the excitation gaps for 
the one doublet modes reaches zero as can be seen in 
Fig. 2(b). The situation is similar to the K° condensa- 
tion (with a small fraction of K + condensation) [39, 40] 
because the ub-sr and db-sg modes belong to the color- 
isospin SU(2)c+y doublet, and have Q = +1 and 0, re- 
spectively. 

When the chemical potential is decreased further, 
the gap parameters split into three different values al- 
though they are all still finite. Accordingly, the isospin 
SU(2)c+v symmetry gets broken by \i e and ^3 so that 
the gaps in quasi-quark dispersions all take different val- 
ues (see FIG. 2(b)). We remark that the electron chem- 
ical potential or its density actually serves as an order 
parameter of the CFL/gCFL (insulator/metal) transition 
as claimed in [11]. 

The gCFL phase continues to be the ground state 
down to /1 = 558.4 MeV below which the UQM phase 
is more favorable in terms of the thermodynamic poten- 
tial. When the transition gCFL — > UQM takes place, 
there should be the large re-configuration of the flavor 
contents as can be seen in FIG. 2(c). Thus, this transi- 
tion requires a lot of electro- weak processes which include 
the d production like ug(ur) — > dg(dr) + e + + v e in ad- 
dition to the decay to the gaplcss modes accompanied 
by the electron production db — ► ub + e~ + v e and the s 
quark decay sb + u — > db + u. We should note that there 
still remains an open interesting question how the UQM 
droplets are dynamically formed in the gCFL phase and 
grow against the surface tension. 

We have seen that, as the density is decreased, the CFL 
phase turns into the gCFL phase, and then the gCFL 
phase gets taken over by the UQM phase. Accordingly, 
the number of gapped quasi-quark modes decrease as 9 
(CFL) 7 (gCFL) -> (UQM) at T = as shown in 
the previous work [24]. Next, we will discuss how the 
situation is changed in the T ^ case. 

b. Phase transitions and crossovers for finite temper- 
ature: We show here that the quark matter undergoes 
a sequence of the transitions, CFL(9) — > gCFLg(8) — > 
uSC(6) -► guSC(5) -► 2SC(4) -► g2SC(2) -► UQM(O) 
as T becomes large; the number of the gapped modes in 
each phase is indicated with a parenthesis. 

In FIG. 3(a), we show the T-dependence of the gaps, 
masses, and chemical potentials. The figure shows that 
the Ai first melts and then A2 disappears as T is in- 
creased. The thermally-robustest pairing is the 2SC, 
while second one is the uSC, which is in agreement with 
the result in [29]. This conclusion will be found to be 
also consistent with the Ginzburg-Landau analysis given 
in Sec. HIE, where we extend the previous study [22] 
and see that the large value of the strange quark mass 
actually disfavors the dSC phase. 

In order to study the thermal transition in more detail, 
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FIG. 3: (a) The gap and mass parameters, and the chemical 
potentials as a function of T at /i = 590 MeV so that the 
system at T = stays in the CFL phase (see FIG. 2). (b) The 
nine quasi-quark gaps versus T. Nonzero Q-charge is also 
indicated by ±1. (c) Just an enlargement of the figure (b). 

we have also plotted the T-dcpendence of the gaps in the 
nine quasi-quark spectra in FIG. 3(b). FIG. 3(c) is just 
an enlargement of (b). The three large points on the 
horizontal line indicate T = T crj (77 = 1,2,3) at which 
A, ; vanishes (the same as the large points in FIG. 3(a)). 
When the temperature is increased from T = 0, the first 
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FIG. 4: The phase diagram for the weak coupling (Ao = 80 MeV) (a), and for the intermediate coupling (Ao = 125 MeV) (b). 



qualitative change occurs at T ~ 15.8 MeV where one 
of the db-sg dispersion becomes gapless and the gCFLs 
phase (gCFL in [31]) sets in. Unlike the gCFL phase 
at T = 0, the ub-sr dispersion does not become gapless 
at this point. From the viewpoint of the insulator-to- 
metal transition, however, the sharp phase transition is 
smoothen not only due to the absence of the gapless ub- 
sr mode with Q = 1, but also to the thermally-excited 
on-shell quasi-quarks at finite temperature as the latter 
ingredient is already noticed in [23]. In fact, from the 
FIG. 3(b), we can see that the ub-sr (Q = +1) mode 
is lighter than the db-sg {Q = —1) mode for T ^ 5 MeV 
and accordingly there is a little excess of the quasi-quarks 
with Q = +1 in the system. In order to achieve the 
neutrality, there must be the equal amount of electrons 
so that fi e takes positive finite value as long as T is finite. 
As a result, the CFL-gCFLg (insulator-metal) transition 
becomes a smooth crossover. 

When T is increased further beyond the gCFLg onset, 
Ai disappears at T ~ 18.4 MeV. This is a second order 
phase transition of gCFLg — > uSC. As a consequence, two 
quasi-quark modes become gapless as is seen in FIG. 3(b); 
one is the sg-db mode i.e., the partner of the gapless bd- 
sg mode, while the other is the isospin singlet dg-sb-ur 
mode. 

Next qualitative change occurs when T is increased 
to T ~ 21.2 MeV. At this point, the ub-sr mode with 
Q = +1 becomes gapless and the guSC(5) phase sets in. 
Notice that because the thermally-excited quasi-quarks 
are already present in the system over the range (E u f,. sr ^ 
T), there is no sharp boundary between the uSC(6) and 
guSC(5) phases. 

At a little higher temperature T ~ 21.8 MeV, A2 van- 
ishes and above which the 2SC(4) phase is realized. At 
this transition point, the sr-ub mode having Q = +1 be- 
comes gapless. Through the gap parameter A2, this sr- 
ub mode is paired with its partner, i.e., the gapless ub-sr 



mode. When A2 — * 0, these two modes get unpaired 
to become the bare ub and sr quarks. This guSC(5) — > 
2SC(4) phase transition is of second order. 

As T is increased further, the crossover 2SC(4) — > 
g2SC(2) takes place at T ~ 26.6 MeV, and finally the 
g2SC(2) is taken over by the UQM(0) phase through a 
second order phase transition at T ~ 27.4 MeV. 

We have found that the excitation gaps behave in a 
more complicated manner than the gap parameters as 
functions of T. However, as we explained, there are no 
sharp boundary with thermodynamical singularity be- 
tween CFL(9) and gCFL 8 (8), uSC(6) and guSC(5), and 
2SC(4) and 2SC(2) transitions. For this reason, we did 
not indicate these crossover boundaries in the phase dia- 
gram of FIG. 1. 



B. Phases for weak and intermediate coupling 

Let us next examine the weak and intermediate cou- 
pling cases (A = 80 and A = 125 MeV). In FIG. 4, the 
phase diagrams for the both cases are shown. 



1. Phases for the intermediate coupling 

We first discuss the case of the intermediate coupling 
Ao = 125 MeV, in advance of the weak coupling case. 

FIG. 4(b) shows that the UQM and gCFL phases dis- 
appear at T = and the g2SC (CFL) phase exist in the 
intermediate (high) density regime [24]. We notice that 
the fully gapped 2SC phase exists in a small /x-region be- 
tween the CFL and g2SC phases. Also we notice that the 
transition between the 2SC and g2SC phases is of a first 
order accompanied by a jump in the dynamical strange 
quark mass and other physical quantities. This is in con- 
trast to the usual 2SC/g2SC transition without strange 
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FIG. 5: (a) The gap and mass parameters, and the chemical 
potentials versus p along the (g)2SC solution. The dashed line 
placed on p = 496.4 MeV represents the point above which 
the 2SC state is metastable against the CFL. (b) The four 
quasi-quark gaps as a function of p,. (c) The fraction of 
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pairing (r, g) sector as a function of momentum p in the g2SC 
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quarks. We will later discuss this point in detail. 

The phases for T ^ also differ from the extremely 
weak coupling case: (i) The window for the uSC phase 
is pushed away to higher density side and is confined 
in a small region. The cross placed on the dashed line 
represents the point at which the window for the uSC 
phase opens, (ii) The large dot putted on the CFL/2SC 
transition line indicates the critical endpoint; the dashed 
line between the large dot and the cross shows that there 
exists the continuous phase transition of 2SC <-> CFL. In 
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FIG. 10: Solid lines are the occupation numbers for the u and 
d quarks in the g2SC phase (p, T) = (500 MeV, 25 MeV) as a 
function of momentum p. The sharp Fermi surfaces for the u 
and d quarks at T = are also indicated by the dashed lines. 
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UQM phase at T = 25 MeV. The inset is the quantity n 
n t , i.e., the isospin mismatch in the UQM at T = (dashed 
line) and for the unstable UQM at T = 25 MeV (solid line). 



FIG. 8: (a) The gap parameter A3 and the chemical potential 
/Lt e , 8 versus T for A = 125 MeV and fx = 450 MeV. (b) The 
excitation gaps in the four quasi-quark spectra as a function 
of T. (c) The SU(2) co i or neutral isospin density divided by 
the total isoscalar density as a function of T. This quantity 
is an order parameter for the 2SC/g2SC transition at T = 0. 



the following, we shall discuss the detail of the 2SC/g2SC 
transitions. We first study the T — case and then 
investigate the T / case. 

a. 2SC/g2SC transition at T = 0: In FIG. 5(a), we 
explicitly show the solution of the gap equation in the 
2SC sector at T = 0. Just to make the argument clear, 
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we show the 2SC solution for the entire ^t-space, but 
note that the 2SC phase is metastable being energeti- 
cally taken over by the CFL phase in the right hand side 
of the vertical dashed line placed on /j, = 496.4 MeV; we 
should keep in mind that the window for the 2SC phase 
is extremely narrow. Below the critical chemical poten- 
tial fjt = 494.8 MeV, (fi ug - \i dr )j2 = {^ ur - fi dg )/2 = 
He/2 becomes greater than the gap parameter A3, and 
as a consequence, the fully gapped 2SC is taken over 
by the partially ungapped g2SC (See FIG. 5(b)). In 
FIG. 5(c), the ratio of the SU(2) co i or neutral isospin den- 
sity J2 a =r g( n d ~ n u) relative to the isospin scaler density 
is shown as a function of This quantity serves as an 
order parameter with which the g2SC phase can be dis- 
tinguished from the 2SC phase; in the fully gapped 2SC 
phase, this quantity must be zero because the equal num- 
ber of u and d quarks are required for the pairing, while in 
the g2SC phase, the quantity becomes nonzero because of 
the accumulation of isospin charge in the blocking region. 
This is clearly seen in the plot of the occupation num- 
bers in the g2SC phase (See FIG. 6). From all the three 
figures (FIG. 5) in addition to the analyses of the thermo- 
dynamic potential, we can conclude that the 2SC/g2SC 
transition is of first order unlike the usual case without 
the chiral dynamics, i.e., the (ss) condensation. 

The first order phase transition between these two 
phases is brought about actually by the competition be- 
tween the chiral (ss) condensation and the quark pairing 
dynamics in the u-d sector as was first recognized in [24] . 
We can see a small jump in the strange quark mass at 
the transition point in FIG. 5(a). The strange quark 
mass reaches almost the vacuum value on the g2SC side, 
while it varies on the 2SC side and becomes smaller as 
the density goes high. Accordingly the strange quarks 
are absent in the g2SC side, while a small number of 
s quarks are present in the 2SC side and it grows with 
increasing density. This is naturally understood as fol- 
lows: In the fully gapped 2SC phase, the u(rg) and d{rg) 
quark densities should be equal for accommodating the 
A3-condensate. Because of the positive electric charge 
1/3 = 2/3 — 1/3 coming from this SU(2) co \ m sector, the 
system needs strange quarks, electrons and dfe-quarks. 
Thus, the high density condition [M s <C n] favors the 
2SC realization because plenty of strange quarks provid- 
ing negative charges can exist in the system as can be 
confirmed by the plot of strange quark density in FIG. 7. 
But this condition is in turn disfavored with respect to 
the chiral condensation energy. On the other hand, the 
g2SC phase can realize the electric neutrality with a less 
number of strange quarks because of a d quark excess to 
u quarks in the SU{2) co \ m sector. Accordingly the g2SC 
phase is not so much disfavored by the chiral (ss) con- 
densation (M s <; /i) in comparison with the 2SC phase 
[60]. In short, the g2SC phase can co-exist with the chi- 
ral symmetry broken phase easier than the fully gapped 
2SC phase can do. The competition of these two ingredi- 
ents makes the 2SC/g2SC transition be of first order at 



T = 0. 

b. g2SC/2SC/g2SC crossovers for T ^ 0: We now 
investigate how the situation changes for the 2SC/g2SC 
transition at T / 0. The phase diagram displayed 
in FIG. 4(b) shows that, in a relatively low density 
regime (430 MeV £ \x <> 480 MeV), the g2SC/2SC and 
2SC/g2SC transitions take place successively when the 
quark matter is heated. This kind of exotic situation is 
already found in [41] where the systematic study of the 
coupling strength dependence of the 2SC/g2SC transi- 
tion is done within the two-flavor NJL model. In our 
case, the strange quark mass reaches almost its vacuum 
value in this region so that the pairing dynamics be- 
comes similar to those in the previous two-flavor mod- 
els. In FIG. 8, we show how the physical quantities be- 
have through these phase transitions for fi = 450 MeV. 
FIG. 8(a) shows the gap parameters and the chemical 
potentials as a function of T. fi e /2 and A3 cross each 
other twice with increasing T so that the system under- 
goes the 2SC/g2SC and g2SC/2SC transitions before the 
pairing is overwhelmed by the UQM. FIG. 8(b) shows the 
gaps in the quasi-quark spectra. We note, however, that 
the excitation gap in the dg-ur (dr-ug) quasi-quark spec- 
trum is always smaller than T so that these quasi-quarks 
are thermally-excited irrespective of whether the system 
is in the 2SC phase or the g2SC phase. These quasi- 
quarks smear the 2SC/g2SC and g2SC/2SC transitions. 
Accordingly, no thermodynamic singularity is associated 
with these transitions except for the final g2SC-UQM 
(A3 — > 0) transition. This can be understood also in 
terms of the order parameter for the 2SC/g2SC transi- 
tion at T = 0. FIG. 8(c) shows the SU(2) cok)r neutral 
isospin density ^2 a=r g {rij — n°) normalized by isoscalar 
density as a function of T. From the figure, we can see 
that the g2SC/2SC transitions are crossovers; the order 
parameter always takes nonzero value and the system 
is always isospin charged due to the thermally-excited 
quasi-quarks even for [i e j2 < A3 region. 



2. Phases for the weak coupling 

Now we turn to the weak coupling case (Ao = 
80 MeV). The phase diagram is shown in FIG. 4(a). In 
this case, the diquark coupling is not stronger enough to 
exclude the UQM phase from T = 0, while the g'CFL 
phase is overwhelmed by the UQM with a large (ss) and 
thus is washed out as was first noted in the previous work 
[24]. 

One of the most striking features in this case is the 
appearance of the g2SC phase for T ^ even when the 
UQM phase is realized at T = 0. This somewhat strange 
aspect of the color-neutral g2SC phase was first recog- 
nized in [41] and has been further confirmed in [23], al- 
though it seems that physics for this phenomenon has 
not yet been clarified enough. In FIG. 9, we show the 
gap parameter A3 and the chemical potential (/i e /2,/is) 
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as a function of T in this situation (/i = 500 McV) . At 
T = 0, the system is in the UQM phase, and the A 3 - 
condcnsate appears at T ~ 16 MeV and then the 2SC 
turns again into the UQM phase at T ~ 36.5 MeV. fi e /2 
is always larger than A3 so that the system is in the 
"g2SC" phase accompanied by only two gapped quasi- 
quarks, i.e., the ur-dg and ug-dr modes. However, it 
should be again noted here that the two phases (2SC 
and g2SC) are not thermo dynamically distinguishable for 
T 7^ as explained in the previous section. To under- 
stand the appearance of the 2SC phase for T ^ 0, we 
show the occupation number of the u and d quarks for 
various situations in FIG. 10. There are sharp Fermi sur- 
faces at T = which are indicated by dashed lines, i.e., 
one for d quarks at ppd = fx + £te /3 and the other 

for u quarks at pp u = f/,— 2/ii T_0 ' ) /3. The dashed curves 
show the u and d quark occupancies at T = 25 MeV in 
the UQM phase which is unstable to the formation of 
u-d pairs. Quarks are thermally-excited so that the dis- 
tribution is somewhat smeared. We note that because 
of this thermal effect, the u and d quarks are under a 
better kinematical matching than the situation at T = 0; 
from the inset of the figure, we can see that the isospin 
mismatch (ud — u u ) at the averaged Fermi momentum 
p = (pFd +Pf*u)/2 is reduced in the T = 25 MeV case 
from the unity, i.e., the value in the T = case. As a re- 
sult, it is easier for the quark matter in the UQM phase to 
form the diquark condensate in the isospin singlet chan- 
nel at finite temperature than at T = 0. The solid lines 
show the occupation numbers after this reconfiguration 
takes place (the g2SC phase). We remark that unlike the 
occupation numbers in the g2SC phase at T — drawn 
in FIG. 6, there are no singular points in those for the 
thermally-smeared g2SC phase for T^0. 



C. Phases for (extremely) strong coupling 

We next discuss the strong and extremely strong cou- 
pling cases, A = 160 MeV and 200 MeV, respectively. 
The phase diagrams are depicted in FIG. 11. The 
phase boundaries move in the directions indicated by the 
bold arrows when the diquark coupling is increased from 
A = 160 MeV to A = 200 MeV. One can see that both 
the phase structures in these cases are qualitatively iden- 
tical and are much simpler than those for the weak and 
intermediate couplings. The premature gapless phases 
disappear from the low temperature regime, and only 
the two major pairing phases remain in the phase dia- 
gram, i.e., the CFL phase and the 2SC phase. In each 
case, we still have a small window for the uSC realiza- 
tion as a precursory phase of the CFL — > 2SC transition 
in the high //-region of the phase diagram. However, 
one should notice that the critical temperature for the 
2SC/UQM transition has a maximum in this chemical 
potential region, which might imply that this shrinkage 
of the uSC window is due to a cutoff artifact. It should 
be also noted that when the diquark coupling is large, 



1 



1 


1 1 


1 1 1 1 

UQM 








X 


\ 






\ 2SC 


/ / CFL 

' / A =160[MeV] 


1 


w 

III ll 


/ A fl = 200[MeV] 
III 



160 
140 
120 
100 
80 
60 
40 
20 


300 350 400 450 500 550 600 650 
H[MeV] 

FIG. 11: The phase diagrams for the strong coupling (Ao = 
160 MeV) and for the extremely strong coupling ( Ao = 
200 MeV). Phase boundaries shift in the directions indicated 
by the bold arrows when the diquark coupling is increased 
from A = 160 MeV to A = 200 MeV. 



the 2SC phase encroaches upon the domain for the %SB 
phase. In addition, the 2SC/UQM transition temper- 
ature becomes larger with increasing diquark coupling, 
and at some critical coupling (Ao = A c ), it reaches the 
tricritical point of the chiral phase transition: At this 
coupling, the chiral tricritical point acquires the nature of 
the doubly critical point where the two second order crit- 
ical lines merge; one is for the 2SC/UQM transition and 
the other is for the UQM/xSB transition. Because this 
point has both the tricritical and doubly critical natures, 
we call it the "TDCP" point. When the diquark coupling 
is increased beyond this critical coupling A c , the TDCP 
point shifts to higher temperature because the 2SC/xSB 
transition is always of first order [29] . Our calculation in- 
dicates the critical coupling A c where the tricritical point 
obtains a DCP nature lies between Ao = 160 MeV and 
A = 200 MeV. 

Finally, it is worth mentioning that the phases in the 
strong coupling cases discussed above are all free from the 
instability problem associated with the imaginary Meiss- 
ner masses [61]. 



D. Baryon density ps at T = 

We have investigated the phase structures of the quark 
matter in the (fi, T)-plane for several diquark couplings. 
It is, however, sometimes more convenient to describe 
the system in the (ps,T)-plane to have a physics intu- 
ition into the system. Also it is necessary to have the 
equation of state as a function of Pb to clarify the inner 
structure of compact stars [30] . Here one should notice 
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FIG. 12: The baryon number density at T = in the unit of 
po = 0.17fm -3 , i.e., the normal nuclear density. 1, 2, and 3 in- 
dicated for each line correspond to the following three cases; 
the extremely weak coupling (Ao = 25 MeV), the interme- 
diate coupling (Ao = 125 MeV), and the extremely strong 
coupling (Ao = 200 MeV), respectively. The cross placed on 
the CFL line indicates the transition point at which the CFL 
continuously turns into the gCFL phase. 



that when multiple phases coexist at a temperature with 
the same chemical potentials, each phase may have a dif- 
ferent baryon density from each other; these phases are 
actually realized in a mixed phase where the phases with 
different baryon densities coexist. 

Here, we ask how dense each phase at T = is. In 
FIG. 12, we show the baryon density pb versus chem- 
ical potential fi at T = in the following three cases; 
(1) the extremely weak coupling Ao = 25 MeV, (2) the 
intermediate coupling Ao = 125 MeV and (3) the ex- 
tremely strong coupling Ao = 200 MeV. In any case, as 
the chemical potential is lowered from the highest value, 
the baryon density decreases and shows several jumps as- 
sociated with first order transitions. We start with the 
case (1). In this case, the ground state at highest den- 
sity, about 15 times the normal nuclear matter density, is 
in the CFL phase. When the density is decreased down 
to the point indicated by the cross, the CFL continu- 
ously turns into the gCFL phase with a steeper gradient 
% = ^(jp) T = k t(p), i.e., a larger compressibility or 
equivalently, a larger number fluctuation (susceptibility 
xj p ( r )) attributed to gapless modes. If the density is de- 
creased further, the first order transition from the gCFL 
to the UQM takes place; accordingly the baryon den- 
sity drops down to almost ten times the nuclear density 
through a mixed phase of the two phases. The UQM 
phase continues down to ~ 3 times the nuclear density, 
and then gets taken over by the chiral symmetry bro- 
ken phase. The situation becomes slightly complicated 



in the case (2). The densest phase is again the CFL. As 
the density is decreased, the CFL phase persists down 
to almost ten times the nuclear density. Then the quark 
system undergoes a first order phase transition to the 
2SC phase via the dynamical nucleation and growth of 
the low density 2SC droplets in the CFL phase. This 
nucleation also needs a lot of weak processes. The den- 
sity window accommodating the 2SC phase is extremely 
small and the 2SC phase gets immediately taken over by 
the g2SC phase. The g2SC phase finally changes into the 
%SB phase. Similarly in the case (3), the baryon density 
drops twice according to the CFL/2SC and 2SC/xSB 
transitions. It is interesting that the lowest baryon den- 
sity of the CFL phase decreases with increasing diquark 
coupling while the density at the same \i is enhanced 
by the factor A^/fi 2 . It should be noted, however, that 
the p-fj, relation will be affected by the fluctuation effect 
so that the density becomes even larger due to bosonic 
degrees of freedom in the strong coupling cases [43] . 



E. The Ginzburg-Landau analysis 

In this section, we make a systematic Ginzburg-Landau 
analysis in order to understand why the dSC phase does 
not appear in our calculation of the phase diagram. Our 
analysis is an extension of the previous work [22]; we 
expand the Ginzburg-Landau coefficients up to quartic 
order in M s which were not taken into account in [22] 
but turn out to play an important role for the relatively 
strong coupling (low density) regime. We also notice that 
this extension makes it possible to provide a unified and 
systematic description of the thermal pairing/unpairing 
transitions obtained in the NJL [29] and the Ginzburg- 
Landau analysis [22]. 

If the phase transition at finite temperature is of sec- 
ond order where the condensate A,, vanishes at the criti- 
cal point, we can expand the effective potential in terms 
of the gap parameters near the critical temperature T c . 
In this section, we focus on the regime where the chi- 
ral symmetry is restored, and assume that the thermal 
pairing/unpairing transitions are not so much affected by 
the chiral dynamics there. Thus, we shall switch off the 
scalar coupling G s in our NJL lagrangian, Eq. (1), and 
treat the strange quark mass M s as a constant parame- 
ter. This is nothing but a model used in [23, 44], which 
will be referred as the diquark NJL model hereafter in or- 
der to distinguish it from the present NJL model analysis 
with the scalar coupling, i.e., Eq. (1). Treating M s as a 
parameter can be justified by the fact that M s does not 
so much vary in the vicinity of T c (See FIG. 3(a)). Ei- 
ther from the diquark NJL model or from the QCD-likc 
theory with the Cornwall- Jackiw-Tomboulis potential up 
to 2PI graphs [22], we can derive the Ginzburg-Landau 
potential. We first emphasize that each of the Ginzburg- 
Landau coefficients becomes a function of M s under the 
neutrality constraints. We evaluate these functions by 
the Taylor expansion in M s . Our final task is to cal- 
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culate the splittings in the critical temperatures up to 
quartic order in M s . Here, we shall only give our result 
for the Ginzburg-Landau potential leaving the detail of 
the calculation to the Appendix A, because it is some- 



what involved although straightforward. After solving 
the neutrality constraints, the Ginzburg-Landau poten- 
tial is found to be 
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where N[fi] = h 2 /2tt 2 is the density of state and T c o 
is the critical temperature for the symmetric (M s = 0) 
quark matter. Some remarks arc in order here, (i) The 
coefficients for the A 4 terms are not expanded in M s . 
This is because those effects on the splittings of T c o are 
small with a suppression factor (T c o//Lt) 4 in comparison 
with the contribution from M s -dependent terms in the 
coefficients of A 2 terms, (ii) Up to quadratic order in 
M s , this exactly coincides with the form obtained in the 
previous study [22] as it should be [62]. (iii) The M 4 - 
correction to the Ginzburg-Landau coefficient is not so 
small because it is not a simple expansion in M s //i; the 
coefficient of M 4 is enhanced by factor /i 2 /T 2 . 

We now calculate the splittings of critical temperature 
T c o — * {Tci, T C 2, Tcz). We first solve the gap equation 
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in A 2 , A 2 , A 2 and define the temperatures (T ds , T su , T ud ) 
as the solutions of A 2 (T ds ) = 0, A^(T SU ) = 0, A§(T ud ) = 
0. Up to the quartic order in M s , we obtain 
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Here we can also see that the M 4 -correction could give a 
significant contribution comparable to the M 2 term be- 
cause of an enhancement factor /i 2 /T 2 log(/i/T c o). Wc 
can see that T ud is largest of the three temperatures, 
which simply means that the 2SC phase is robustest 
against the thermal disturbance. We can evaluate the 
critical temperature T c3 for the 2SC — > UQM transition 
by putting Ai = A 2 — > in Eq. (24) and solving = 



in T. The result is 
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In contrast to the fact that the 2SC is always the hottest 
pairing phase, the second pairing phase next to the 2SC 
depends on the value of the strange quark mass. We have 
to consider the following two cases; 
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In case [A] , A 2 first becomes zero at T su = T c2 and then 
Ai vanishes at T c \(> T c2 ) when the quark matter is 
heated from the CFL phase; that is 
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This is actually the case which is studied in [22]; the 
quark matter undergoes a hierarchical unlocking CFL — > 
dSC — ► 2SC — > UQM. We note, however, that the quartic 
terms are derived for the first time in this study and these 
will turn out to be a crucial for a unified picture of the 
thermal unpairing phase transitions. 

Let us now examine the case [B] with a relatively large 
strange quark mass. In this case, Ai first melts at T c \ = 
Tds and after that, A 2 disappears at T c2 (> T c i); that is 
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FIG. 13: The plots of T CI) evaluated by the Ginzburg-Landau 
analysis with typical parameter set Ao = 25 MeV and /j, — 
500 MeV. When the coupling strength is increased, the criti- 
cal DCP stress moves to higher value according to Eq. (34), 
while the ratio of T c o — T C 3 : T c q — Tdcp = 1:7 would not be 
affected within the approximation we are working. 



In this case, the quark matter undergoes another hierar- 
chical unlocking CFL -> uSC 2SC -> UQM. 

We have seen that the dSC phase is the second hottest 
pairing phase in the case [A], while the uSC phase takes 
over the dSC phase in the case [B] . One can see that the 
former (latter) situation will be realized at high (low) 
density. In fact, the case [A] includes the high density 
situation with the scale hierarchy 
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while the case [B] includes the low density regime with 
the scale hierarchy 
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Thus the case [B] contains an interesting regime with 
A ~ where the window for the gaplcss phases 

open at T = 0. These two physically distinct regimes are 
separated by the doubly critical point (DCP) which was 
first discussed in the numerical analysis of the diquark 
NJL model [23]. We can analytically derive the critical 
stress for the DCP point in our current framework as 
follows; 
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Just at the density corresponding to this DCP point, we 
have a simultaneous melting of Ai and A2, and hence a 
direct transition CFL — > 2SC takes place with increasing 
T [23]. We can express critical stress for the DCP in 
terms of the CFL gap A at T = and M s = with the 
aid of the weak coupling universal relation 
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Thus, the DCP point can be parameterized by the cou- 
pling strength Aq as 
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Also, the ratio of T c q to the CFL — > 2SC transition tem- 
perature at the DCP point which we denote by Tdcp = 
T c i = T C 2 can be evaluated as 
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is universal being independent of the coupling choice Ao. 

FIG. 13 shows the phase diagram calculated with 
the Ginzburg-Landau potential for Ao = 25 MeV and 
(i = 500 MeV, where the critical temperatures (T c „) ver- 
sus the stress energy (M 2 /fi) is depicted by changing 
M s with fj, fixed. The DCP is located at (Mf /m)dcp = 
26.6 MeV < n = 500 MeV and T c DCP /r c0 = 0.90. These 
values agree well with those obtained in the numerical 
analysis of the diquark NJL model where (M 2 /V)dcp ~ 
30 MeV and T r DCP /T c0 - 0.85 [23]. If we increase the 
coupling to A = (40, 100) MeV, then the DCP shifts 
to (M s 7/i) DC P = (58.5, 248.4) MeV, and T C DCP /T c0 = 
(0.81,0.44) according to Eqs. (32) and (35). We can 
see that the agreement with the diquark NJL result 
[23] becomes worse in the stronger coupling. This is 
because MP CP shifts to higher value as Ao increases 
so that the Ginzburg-Landau approach becomes worse 
due to lack of the scale hierarchy M s <C /i. In TA- 
BLE II, we have summarized the parameter regime where 
the Ginzburg-Landau approach and the diquark NJL 
model are valid. It is worth stressing that according to 
the Ginzburg-Landau evaluation of the DCP given by 
Eq. (32), (M s 2 /^)dcp is located around - A (A //z), 
which is lower with a factor A c o//i than A c o, i.e., the 
value corresponding to the interesting low density regime 
Mil [i ~ A ; the dSC phase is realized at higher density 
than the gCFL/CFL transition density. 

In the following, we discuss where in the (//, T)-plane 
the DCP is located. We first discuss the possibility that it 
is in the moderate density regime using the diquark NJL 
model. Since the full thermodynamic potential should 
take the form like Q((M,M S ) = T C 4 /(M S /^, T c0 /m) and 
T c o//x is a slowly varying function of /i, changing M s and 
changing fj. are independent of each other in principle. If 
we fix M s to some value, we can obtain the phase diagram 
in the (/j,, T)-plane like FIG. 2 of [44] by determining the 
phase by changing (/a, T); this is contrasted to FIG. 13 
where M s is changed to control the density. In this case, 
we can find the unique DCP in (/^,T)-plane, which we 
denote by (/xdcpi^dcp) = (Mj^)dcp- We note that the 
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Regularization scheme 


Theoretical parameters 


Parameter regime of validity 


Ginzburg-Landau theory 


mass counter term 


T c0 , M s , /j, 


T ~ T c0 , n » M a ,T c0 


high density 


Diquark NJL model 


momentum cutoff 


A, Gd, M 3 , ft 


A ^ fi ^ h x sb 


moderate density 



TABLE II: The comparison of the Ginzburg-Landau theory and the diquark NJL model, i.e., the NJL model without the 
scalar (gq)-coupling [11, 23, 44]. In the derivation of the Ginzburg-Landau potential, we adopt the regularization via the mass 
counter term (Thouless criterion), which insures that all the dimensionful mass scales vanish (correlation length diverges) at 
T c o when M a = 0. This regularization is model-independent in the sense that, after the prescription, the effective potential have 
neither A-dependence nor the explicit Gjj-dependence. The Ginzburg-Landau approach is valid for the weak-coupling regime 
of QCD at high density 3> M s ,T c o), while the diquark NJL model is expected to describe the non-perturbavative regime of 
(A <; n /x x sb), where /x x sb is the chemical potential for the chiral symmetry restoration. 



DCP (/x,T)dcp shifts in (/i,T)-plane if M s is changed. 
We have plotted the DCP (/z,T)dcp for several choice 
of M 8 in FIG. 14(a). Although we have used the same 
diquark NJL model as [44] , we did not adopted the chem- 
ical potential shift approximation as is made there, but 
performed the exact treatment of M s since the condition 
AI S <C [i does not hold well. As can be seen from the 
figure, as M s is increased, the DCP shifts to larger /x 
pushing away the dSC phase to higher density. We did 
not find the DCP for M s ^ 187 MeV, where we find only 
the uSC in the phase diagram. This may be a cutoff 
artifact because /idcp increases with increasing M s and 
approaches A = 800 MeV. 

We have seen that the location of the DCP in the 
(/i, T)-plane has one-to-one correspondence to the value 
of M s . Conversely, if the DCP (/i,T)dcp is found for 
some fixed strange quark mass M s , then this M s can be 
viewed as the doubly critical strange quark mass M^"~ p 
for fixed chemical potential (i = /^dcp; this follows from 
the observation that the dSC (uSC) phase will be real- 
ized for M s < Af s DCP (M s > Af s DCP ) as can be expected 
from FIG. 13. In this sense, we have plotted the doubly 
critical strange quark mass M S DCP as a function of /j in 
FIG. 14(b). For Af s -region above (below) Mf CP (ri) line, 
the uSC (dSC) phase would be realized. Just for com- 
parison, we have also shown the weak coupling Ginzburg- 
Landau evaluation of the doubly critical strange quark 
mass, which can be obtained via putting (/i, T c o)-relation 
into Eq. (32) (sec solid line indicated by Af s DCP (/i, T c0 )). 
We can see the sizable deviation between these two lines, 
which is simply attributed to M s <^fj,. However, we stress 
that the qualitative behaviours are the same; if one goes 
higher density, the doubly critical strange quark mass 
shifts to larger value so that the dSC phase becomes ro- 
buster against the strange quark mass. In order to un- 
derstand why the dSC does not appear in our model with 
the scalar coupling G s and (gg)-condensates, we have also 
plotted the dynamically determined strange quark mass 
in the xSB/UQM sector (see the solid line indicated by 

M^ SB ). The dynamical mass is evaluated on the 
/i-T c0 line in FIG. 14(a) so that it should be regarded 
as the lower limit of the dynamical strange quark mass 
in the superconducting phase. We can see that the line 
of M s xSB/UQM is well above the line for A/ S DCP (^) and 
the two lines never intersect for ll A, which provides 
one probable reason for the abscense of the dSC in the 



present NJL model with dynamical chiral condensates. 

Let us finally discuss the possibility that the DCP is 
located in the weak coupling (/i, T)-regimc of QCD. In 
QCD, it should be unique because the strange quark mass 
M s is a decreasing function of ii approaching its current 
value m s ~ 100 MeV in the high density limit, while the 
doubly critical strange quark mass is an increasing 

function of il. If we assume that the crossing (DCP) point 
/j, = ll* at which A/ S DCP (^*) = M s (n*) is located where 
M° CP (/i*) <C [i*, then we can estimate [i* by the weak 
coupling Ginzburg-Landau result, Eq. (32). By using the 
universal relation, Eq. (33), and the weak coupling per- 

turbativc formula for gap, Ao(/i) ~ [ig e [45 ; 46] 
as well, we have the formula M® CP ~ Ao(/x)/\/g; this 
is actually a slowly increasing function of [i when g is 
identified with the running coupling constant g(pL) vary- 
ing with /i. According to the Schwinger-Dyson analy- 
sis of the gap [46-48], A (/i) / 'y/^tp) reaches - 100 MeV 
around ll ~ 10 10 MeV; this indicates n* would be located 
at even higher chemical potential since m s ~ 100 MeV is 
the lower limit of the mass function M s (fi). It should be 
also noted that our starting assumption Af,p cp (ll* ) <C /i* 
can be justified because the weak coupling condition 
A < A (/i*)/ v /.g(^*) -C pL* is satisfied. We can con- 
clude that the dSC phase is realized in the extremely 
high density regime of fi <; ll* of QCD. 



IV. SUMMARY AND OUTLOOK 

In summary, we have investigated the QCD phase di- 
agram with a special attention to the interplay between 
the chiral and diquark dynamics for a wide region of the 
diquark coupling strength. Our results for the two lim- 
iting cases, the intermediate and strong coupling, quali- 
tatively agree with those obtained in the recent analyses 
[29, 30]. Our central results can be summarized as fol- 
lows, (i) As the diquark coupling is increased, the phase 
diagram gets gradually dominated by the fully gapped 
phases, while premature gapless phases get excluded out 
as is noted for T = in the previous work [24] . In partic- 
ular, the phase diagram in the strong diquark coupling 
(Ao 160 MeV) docs not include gapless phases and 
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FIG. 14: (a) The DCP point (ai,T)dcp for various value of M s ; from left to right, M s = 100 MeV, 110 MeV, • • • , 180 MeV and 
M a = 186 MeV. The quark-quark coupling G d is fixed so that A = 40 MeV at fi = 500 MeV as in [44]. The solid line is 
T c o(/i), i.e., the critical temperature for M s = as a function of fx. (b) The doubly critical strange quark mass M° CP as a 
function of jj,. The solid line indicated by M DCP (/i, T c o(/x)) represents the weak coupling Ginzburg-Landau evaluation of M, DC 
by putting the numerical p-T c o relation depicted in (a) into Eq. (32) 

determined strange quark mass in the xSB/UQM sector of our model with the scalar (gg)-coupling 



The solid line indicated by M * S U Q M j s ^ ne dynamically 



hence the system is automatically free from the instabil- 
ity problem associated with imaginary Meissner masses, 
(ii) The second hottest pairing phase next to the 2SC 
phase depends on the quark density (and also on the 
diquark coupling strength). In particular, the dSC is 
not realized for all the parameter regions for which we 
have performed the calculation. This can be nicely un- 
derstood by the Ginzburg-Landau analysis incorporating 
the quartic-order (Mf) effects to all the pair susceptibili- 
ties. In this analysis, we derived the analytical expression 
for the doubly critical point (DCP) [23] at which the two 
second order critical lines intersect; one for the 2SC/dSC 
transition in the high density side and the other for the 
2SC/uSC transition in the low density side. On the basis 
of this analysis, we conclude that the window for the dSC 
phase at high density (T c0 , A ) S> M s tends to shrink as 
the density is lowered, (iii) We have demonstrated how 
and why non-trivial first order transition could be caused 
by the competition between the diquark and chiral con- 
densations, (iv) We have studied how the 2SC/g2SC 
and the CFL/gCFLg transition is smeared by the ther- 
mal disturbance by regarding the electron density and 
color-neutral isospin density as the order parameters. 

In this work, we ignored a Kobayashi-Maskawa-'t 
Hooft six quark interaction [38, 49, 50], whose effect on 
the pairing is in part taken into account in [29]. This 
term, however, also brings a non-trivial interplay be- 
tween the chiral and diquark dynamics through a mixed 
contraction like, |A3| 2 M S , for example; these terms in- 
deed will be induced by the procedure of the mean field 
approximation. Also we have neglected a possible K° 
condensation in the CFL phase which may be relevant 
to the instability problem associated with the gapless 



phases [51, 52]. Examining this would require a further 
improvement of our understanding on how the meson 
spectroscopy on the CFL will be modified by the neu- 
trality constraints [53]. Extending current work so as 
to take these into account remains indeed an interesting 
future problem. 

A number of papers have been devoted to the reso- 
lution of the imaginary Meissner masses in the gapless 
phases [13, 14]. These include; (i) exploring a possibil- 
ity of the mixed phase [15] along the line suggested ear- 
lier [54-56], (ii) the baryon current generation [16], (iii) 
examining the crystalline pairing phases [57] taking the 
neutrality constraints into account [17, 18], (iv) investi- 
gating a possible secondary gap formation [19], (v) the 
gluon condensation in the g2SC [20] , and (vi) an inhomo- 
gencous (p-wave) Kaon condensate [21]. Examining all 
these possibilities or searching for a more stable ground 
state as well as exploring a nature of the gluonic instabil- 
ity in the gapless phases is now one of central problems 
in QCD. Apart from this instability problem, we have 
shown that the chiral dynamics in the strange quark sec- 
tor may simply wash out all the unstable regime of the 
gapless phases if the strength of the diquark coupling is 
large enough. It should be noted, however, that the /i- 
window for g2SC phase still remains at T = in the phase 
diagram even for the diquark coupling Gd/G s = 0.5, i.e., 
the value extracted from one gluon exchange vertex. De- 
termining the strength of the in-medium diquark cou- 
pling from the phenomenological studies including quark 
models applied to the exotic baryon spectroscopy as well 
as the lattice QCD simulation [58] will be needed for a 
realistic description of the pairing dynamics at low den- 
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sity. 
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where E{ = ^/p 2 + Mf and 
n (Mi,fx ai ,T) = -2T J log (l 

dp log(l + e-^+^)/ T ) 



-(Ei-» ai )/T 



-2T 



(2tt) s 



Because H36 is block-diagonalized as W36 = Hf^ 8 © 

E!S ) 9) ' (sr ' U&) ' ( " S ' dr) [(K /5 ) © (-<)], we can evaluate 
the quadratic and quartic terms for each sector separately 
as below. 



a. (TL^ 13 ) © {-7-Q 13 ) sector: For the symmetry rea- 
son, we only have the term which is even order in A. Thus 
we have to evaluate the quadratic and quartic terms in 
Eq. (Al). The quadratic term can be regarded as the 
static susceptibility. For example, we define x {db ' sa) for 
the (a, pi) — (db, sg) sector by 



APPENDIX A: EVALUATION OF 
GINZBURG-LANDAU COEFFICIENTS 

In this section, we present the detail for a derivation of 
Eq. (24). We work with the NJL model for a while, but 
the results will be found to be model-independent by a 
subtraction of the Thouless criterion: By this subtraction 
scheme, the model dependence is reduced only to the 
critical temperature T c q. 

We first expand Eq. (5) in terms of A at T ~ T c q. Up 
to the quartic order in A, we obtain 
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where we have defined 



dp 



(27T) 



tr 



1 



+ (iuj n <-> —iuJn) 



Alx^t S op(0,0^,T,^fi 3: fi a ,M s ). 



(A3) 

We can analytically evaluate Xuioop ^ expanding it in 
t= ( M s '/fiaP^Hap) with fi a p = (fi a + up) /2 and 8^ a p = 
(fx a — np)/2 being referred as the averaged and relative 
chemical potentials, respectively. Up to quartic order 
in t, we obtain the following expression by performing 
all the Matsubara summations and the energy integrals 
except for the term suffering from the UV divergence, 



-^l-loop 
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(A4) 



H° 6 (p;M,/x) = W 36 |a^o, V 36 (A) = H 



36 



n/0 

H 36- 



We omitted the term for the mean-field potential of M s 
in Eq. (Al) because we treat it as a constant and as if an 
external parameter near T c . The zero-th order term in- 
volves the free quark contribution to the thermodynamic 
potential and the vacuum fluctuation which we shall sim- 
ply subtract. 



- T E/ 7^ trL °g M-^] 



To obtain this, we have ignored the contribution from 
quasi-antiquark poles. Also we have defined the density 
of state by N[fjt] = ^ 2 /2n 2 . 

The term quartic in Ai can also be derived within the 
same approximation as done in the quadratic term: 



N\jl 



dbsg\ 



1 - 



3M 2 \ 7C(3)Af 



4/W ^ 2 T 2 



(A5) 



We note that order t correction corresponding to the sec- 
ond term in the above expression gives rise to an order 
(A2) correction to the critical temperature T cv . However 

it can be shown that this correction to T cr) is suppressed 
with a factor 0(T 4 o //i 4 ) in comparison with that from 
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the t 2 terms in the quadratic term in Eq. (A4). There- 
fore we shall ignore the second term in Eq. (A5) in the 
discussion below. 

We have demonstrated the expansion in the case of 
(db, sg) sector, but the Ginzburg-Landau coefficients in 
other two sectors, (sr, ub) and (ug, dr), can be derived in 
the completely same manner. 

b. H^ s sector: Now we take care of the three- flavor 
mixed sector Hi2 S ■ The terms quadratic in A can be 
shown to have the completely same form as that in the 



(7^4^)9) (— T~Q P ) sector. For instance, the term quadratic 
in A\ becomes 



,a/3 



A 2(dg,sb) 
"^lAi-loop i 



(A6) 



where X if Sop ^ s defined by Eq. (A4) with the replacement 
(db, sg) — > (dg, sb). On the other hand, the quartic term 
coming from this sector turns out to be of the exotic 
form: 

7C(3) 



(Ai; + Ai + A^ 
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8w 2 t 2 

Ml 21C(3) 
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t^uds 
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A 2 



/z sb )/3. Also 



where we have defined Ji u ds = + Hdg 
M 2 correction can be safely ignored in the analysis below 
for the same reason as given previously in the non-mixed 
two flavor sector. 



Using this, we obtain 



(A 2 + A 2 + A 2 )4AV] 
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This is exactly of the same form as that obtained earlier 
in [8, 59]. 

We now extend this result with the effects of strange 
quark mass and the charge neutrality constraints taken 
into account. Our task here is to calculate the splittings 
of the melting (A r) — > 0) temperature T cl) up to the quar- 
tic order in M s . For this purpose, we need the expansion 
of the Ginzburg-Landau coefficients in terms of variables 
t = (M 2 / jjL, ju e , /i3, /ig). In order to derive T cv up to the 
quartic order in M s , we have to expand the coefficient 
of the quadratic term up to second order in t. Thus we 
start with the following Ginzburg-Landau potential. 



c. Ginzburg-Landau potential: We first investigate 
the ideal case of M s = [i e ,3.8 = 0. Combining all the 



results above and putting M s 
have 
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In this case all the susceptibilities xC\ol p become of 
the same form and we denoted them by xi_i op- 

The 

quadratic term is divergent and needs a renormalization. 
We renormalize it with the use of the gap equation at 
T = T c q (The Thouless criterion for the critical tempera- 
ture for three massless quark matter) , which plays a role 
of the mass counter term in [1] : 
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Subtracting this from the quadratic term, we obtain 
' 4 

Gd 



A 2 



2xi-ioo P (0, 0;n,T) 



AA 2 N[iAjdE± 
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tanhj ) 
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(A9) 



4/Grf in the quadratic terms can again be replaced by 
2xi-ioop(0, 0, /i, T c o) by using the Thouless criterion in 
the symmetric matter Eq. (A8). By doing this, the 
model dependence disappears and only its remnant is 
condensed into the parameter T c q. After this replace- 
ment, we can expand the coefficients of A 2 , A 2 and A 2 
in t = (A/ 2 /^i, /i e , ^ 8 ). For example, the A 2 term can 
be expanded up to second order in t as follows. 
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We have used the identity AT [/i + Sfj] = N[/j]+ hS^i/it 2 + 
S/j, 2 /2tt 2 . Due to the asymmetries between cross-species 
which is caused by M s and {// e , A*3, Ms}; we cannot com- 
pletely eliminate the divergent quadratic term by the 
subtraction of the gap equation at T c0 and we still have 
divergent corrections proportional to 
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dE — tanh 

E \2T c0 



log 
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with 7 being the Euler constant and io c is a UV cutoff 
for the quasi-quark energy. We note, however, that this 
divergence is originated simply in our constant gap pa- 
rameter ansatz and can be made finite once the momen- 
tum dependence of the gap parameter is properly taken 
into account; it can be shown, in the same manner as 



adopted in [22], that we can replace the above divergent 
contribution by 



A 2 log 



2w c e n 
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where A^ here should be regarded as the value of the gap 
energy at the Fermi-momentum p = /i, and g is a running 
gauge coupling constant evaluated at energy scale /i. In 
the following, we use the above formula instead of (in- 
dependent expression. The A| and A 2 , sectors can be 
evaluated in the similar way. After combining all the 
results, we have the full Ginzburg-Landau potential up 
to the quadratic order in (Af 2 /^,/i e ,/i 3 ,yLtg): 
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This is one of the central results from which the analytical with 
expression for the splittings of the critical temperature 

(?co — * T cn ) can be derived (see Sec. HIE). We now PFu = A* — 3"Me, 

impose the charge neutrality constraints by solving p pd _ ^ _|_ 

= 0, — — = 0, — — = 0, (A16) PFs = JQi + lnJt-M*. 
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We here come back to Eq. (A16), and solve these equa- 
tions up to the quadratic order not only in Al 2 /fi, but 



also in A^. Wc obtain 



in qwr.i = T^4/T,) (t) 2 - This ma y cx P lain 
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Substituting these expression to the free potential Vt e + 
^2 ai and omitting terms which is independent of A^, 
we have the following A 4 contribution 
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Also the feedback contribution to the quartic term comes 
from the quadratic A 2 terms in Eq. (A15). Substitut- 
ing the chemical potentials to the quadratic terms in 
Eq. (A15) leads to the following quartic terms. 
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These feedback terms bring about the M 2 (M 4 ) correc- 
tions to the critical temperatures T CI; , but these correc- 
tions are C(T ( f //x 4 ) suppressed to those from M 2 (M 4 ) 
terms in the susceptibilities. Therefore, wc can safely 
ignore these feedback contribution as long as [i ^> T c0 . 
Substituting Eq. (A20) into the quadratic and quartic 
terms of Eq. (A15) and picking the parts which survives 
under (i 3> T c q lead to our final result for the Ginzburg- 
Landau potential, i.e., Eq. (24). On the basis of this po- 
tential, we can argue how the CFL pairing gets dissolved 
when T c o is approached as is done in Sec. Ill E. 

In the present framework, we can also give the analyt- 
ical formula for the following quantity, 
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We have derived this expression starting with the 
Ginzburg-Landau potential. Conversely, if the Ms- 
dependence of this quantity is known, it gives some infor- 
mation about the Ms-dependent terms in the Ginzburg- 
Landau potential. In [23], this quantity is expanded 

- A^il^T o0 = « + K^) 2 + <^) 4 + - , and 
the coefficients, a, b and c, are extracted from the nu- 
merical result of the diquark NJL model; a = 2.52, 
b = 36.2 and c = 1.02 x 10 3 are obtained with a pa- 
rameter choice ([i, A ) = (500, 25) MeV. We find here, 
however, that Eq. (A21) does not take a form of a simple 

expansion in (— ) but rather seems to be an expansion 



APPENDIX B: OFF-DIAGONAL COLOR 
DENSITIES 

Here, we show that the off-diagonal color densities au- 
tomatically vanish for the standard ansatz with the di- 
quark condensate, i.e., Eq. (7). To prove this, we include 
the chemical potentials not only for the diagonal but also 
for off-diagonal color charges as given in Eq. (3). The 
Nambu-Gor'kov Hamiltonian density in this case takes 
the following form after spin-degeneracy removing, 
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where all the matrix elements are 4x4 matrices defined 
by 
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respectively. We define by putting /xi = fi2 = = 
M5 = = M7 = in Eq. (Bl). 
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This Hamiltonian density H® 6 can be further reduced 
to the block-diagonalized form as given in the text after 
some unitary transformation which makes the Nambu- 
Gor'kov doubling explicit in the (ug,dr), (db,sg) and 
(sr, ub) sectors. However this is not necessary in the fol- 
lowing argument, so we proceed further with Eq. (B3). 



23 



We now define the off-diagonal color charge matrices 
in the Nambu-Gor'kov bases as 



dn 36 

dfx a 



(B4) 



with a = 1, 2, 4, 5, 6, 7 in addition to the diagonal charges 

dH 36 



Q 



e,3,8 



e.3.8 
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What we have to show is that the off-diagonal color 
densities automatically vanish on the ground state deter- 
mined by 7^36- To see this, we first define the complete 
sets by the following eigen-value equation: 

H° 36 \p,a,a) = ae a {p)\p,a,a)°. (B6) 

We have 36 eigenvalues which we distinguish by a = 
1, 2, • ■ • 18 and the Nambu-Gor'kov spin a = ±. The off- 
diagonal color densities can be written in the form as in 
Sec. II B, 



Pa = -|Ea 8 =l/(^ tanh (lf) °(P' a i+\Qa\P,a,+)° 
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(B7) 
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The Q charge matrix in the Nambu-Gor'kov bases 



Q = -Q e - Q 3 - \q% 



(B8) 



commutes with the Hamiltonian density Tl 3e , i.e., 
[H 36 ,Q] = 0, so that the quasi-particle states \p, a, ±}° 
can be chosen to be the eigenstates of Q. Also we note 

Qi = -i[Q,Q 2 ], 
Q 2 = +i[Q, QJ, 

Qi = -i[Q,Q 5 ], 
Q 5 = +i[Q, Q 4 }. 



(B9) 



From these commutation relations, we can immediately 
conclude 

Pi = P2 = pi = P5 = 0, (BIO) 

because 

o&a^lQ^a,?) =i°{p,a,a\[Q 2 , Q]\p,a,a) a = 

(Bll) 

This conclusion follows from Q is diagonal in the quasi- 
particle bases, 

V a, a\Q\p, (3, a') = q°6 a p6 aa , , (B12) 

where q a a takes the values (—1, 0, +1). 



In contrast, the disappearances of p% and p? cannot be 
proven in the same manner because of the commutation 
relations: 



i[Q, Q 6 ] = 0, i[Q, Q 7 ] = 0. 
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Thus we try to give a more direct proof here. We first 
re-write Eq. (B7) as 
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Since the trace does not depend on the base, we can 
evaluate the trace with the natural base instead of the 
quasi-particle eigen-spinors. Thus what we have to show 
here turns out to be the proof of 

tr [H° 36 Q 6J ] = 0, tr [H° 36 n° 36 n a 36 Q 6 j] = 0, ■ • ■ . (B15) 

If this infinite series of conditions can be shown to be 
true, then we can conclude 



Pe = P7 = 0. 



(B16) 



Using the explicit form of the matrices 7i. 3e and Q 6 7 , 
we have explicitly checked using the M&them&tica, that 
the equation 



tr [(O 2 -^^] = 



(B17) 



indeed holds for n = 1, 2, • • • , 18. We have not confirmed 
Eq. (B17) for n > 18, but the confirmation up to n < 18 
is adequate for the reason we shall give in the following. 
First, we write the Hamiltonian density 



is 
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ct=l 



^36 = XVe«(p).Rra = £ (fa - P- a )s a (p), (B18) 
cr.Q 



where P aa is the projection operator a, cr) 00 (p, a, a\ 
in the natural base; this operator projects vectors out to 
the eigenspace in which H 36 = cre a (p) holds. Using this 
decomposition, we can re-write the 18 conditions as 



= £ a 8 =i^tr[(P Q -P_ Q )Q 6i7 ], 
= Ea 8 =i4tr[(P Q -P- Q )Q 6 , 7 ], 



(B19) 



18 : = 4 5 tr [(P a - P- a )Q 6 , 7 ] ■ 



If all the eighteen roots {e a } with a = 1, 2, ■ ■ ■ ,18 take 
different values, then the above eighteen conditions sim- 
ply mean 



tr [(P a - P- a )Q 6 ,r] = °' 



(B20) 
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for a = 1, 2, • ■ • , 18. Therefore for the arbitrary integer 
n, tr[(H^ 2 ) 2 " _1 Q 6 7 ] = should hold. In the case that 
the degeneracy is present as £\ — £2 = £3, for example, 
we can prove the following equation in the totally same 
manner as in the above argument, 



(B21) 



a=l,2,3 



and again reach the same conclusion tr[(7i3 2 ) 2 " 1 Q 



0. Thus, we have proven that the following condition 
indeed holds for a — 6 and 7. 



tr 
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Q, 
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Consequently, we reached the fact that all the off- 
diagonal color densities automatically vanish under the 
assumption of the diquark condensate given in Eq. (7). 
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